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Abstract 

Motivated by recent advances in the apphcation of effective field theory techniques 
to hght nuclei we revisit the problem of electron-deuteron scattering in these ap- 
proaches. By sidestepping problems with the description of electron-nucleon scat- 
tering data in effective field theories, we show that the effective theory expansion 
for deuteron physics converges well over a wide range of momentum transfers. The 
resultant description of the physics of the two-nucleon system is good up to virtual 
photon momenta of order 700 MeV. 



PACS nos.: 12.39.Fe, 25.30.Bf, 21.45. +v 



1 Introduction 



Electron scattering from nuclei has a long and rich history. In impulse approximation 
the charge form factor probed in such experiments is the Fourier transform of the nuclear 
charge distribution, and so these measurements have often been regarded as independent 
tests of models of nuclear structure [T1I2] . In particular, the structure of nuclei with A < 10 
can now be calculated ab initio from a given two- (and three-) nucleon interaction 
Calculations of electromagnetic form factors of these nuclei then reveal agreement with 
experimental data that is, in general, very good [U Ej 

Here we focus on the simplest non-trivial nucleus: deuterium. Elastic scattering of 
unpolarized electrons from deuterium results in an 0{a'^) differential cross-section: 
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where 6e is the electron scattering angle in the centre-of-mass frame of the collision, 
= {Pe—PeY = —Q^ is the (negative) virtuality of the (single) photon exchanged between 
the electron and the nucleus, and is the Mott cross-section for electromagnetic 
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scattering from a point particle of charge |e| and mass Met. 

Deuterium is a spin-one nucleus and so has three independent form factors. These are 
usually denoted by Gc, Gq, and Gm- They are related to Breit-frame matrix elements 
of the deuteron electromagnetic current, J^, through: 
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where we have labeled the deuteron states by the projection of the deuteron spin along 
the direction of the three- vector — pe, and r] = Q^/{AMj). When defined in this way 
these charge, quadrupole and magnetic form factors have the normalizations: 
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where Qd = 0.286 fm^ [6j is the deuteron quadrupole moment, and /i^ 
deuteron magnetic moment in units of nuclear magnetons. 
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The experimental quantities A and B can then be computed from theoretical models 
of deuterium, since 
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However, it was not until the development of experiments with polarized deuterium targets 
that it became possible to unambiguously extract both Gc and Gq from electron-deuteron 
scattering data. The tensor-polarization observable, T20, is related to the ratios 
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and so a measurement of T20, together with measurements of A and B allows an ex- 
traction of Gc and Gq, and hence a complete test of our theoretical understanding 
of deuteron structure. Experiments over the last dozen years at Bates |H1 E], Novosi- 
birsk pm ITT] , NIKHEF ^21 ^1 5 Jefferson Laboratory ^H] have measured T20 in 
electron-deuteron scattering, and so facilitated experimental determinations of the full 
set of deuteron structure functions over a kinematic range between Q = and Q = 1.5 
GeV Modern nucleon-nucleon potentials, when combined with models for two- 

body contributions to the deuteron current, do a good job of reproducing this data (see, 
e.g. [IHl Uni EH! ) • ^ thorough status report on the subject of electron-deuteron scat- 
tering we refer to three recent reviews which discuss the subject [T71 EH 122] ■ 

In this paper we wish to address electron-deuteron scattering data in the frame- 
work of effective theories of deuteron dynamics. This approach (for recent reviews see 
Refs. j2nil21|) is based on the use of a chiral expansion for the physics of the two-nucleon 
system. Ultimately it shares many features with the more "traditional" , and very success- 
ful, potential models. However, as first suggested by Weinberg [2312011221, this "nuclear 
effective theory" is based on a systematic chiral and momentum expansion for the the ker- 
nels of processes in the NN system. Thus, for electron-deuteron scattering we expand the 
deuteron current in operators which are ordered according to their chiral dimension, 
viz.: 

00 -1 

^E^^x^c)«, (13) 
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where the operator (9^ contains i — 1 powers of the small parameters p (the momentum 
of the nucleons inside deuterium), m^, and Q. The numbers Cj are, a priori, assumed to 
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be of order 1, and A is tlie scale of cfiiral symmetry breaking: A ~ 47r/^, mp, M. Since 
tlie expectation value of p and the value of are both much smaller than A it follows 
that, provided Q < A, and the cs really are of order one, this expansion should converge 
well. The expansion parameter (p, (5,m^)/A is denoted here by P. 

The operators O^^'' and the coefficients q are constructed according to the well-established 
counting rules and Lagrangian of chiral perturbation theory |2H1- Here we present results 
for Gc and Gq up to order eP^, and results for Gm up to Oi^eP"^). We go beyond the 
recent calculation of Ref. [22], which computed all three form factors only up to O(eP^). 
We also demonstrate that, provided single- nucleon structure effects are correctly included 
in the calculation, the nuclear effective theory is, in fact, much more accurate than the 
results of Ref. [211 might lead one to believe. Indeed, ultimately it describes all of the 
extant experimental data on Gc and Gq out to momentum transfers of order 700 MeV. 

This is done as follows. In Section |21 we sketch the derivation of from the counting 
rules of chiral perturbation theory, and give results for the current at leading order, 0(e), 
next-to-leading order O(eP^), and next-to-next-to-leading order, 0{eP^). In Section [HI 
we will discuss the wave functions used in our calculation, and outline some of the issues 
associated with the desire for consistency between the deuteron current and the deuteron 
wave functions. In Section [H we will present our results for Gc, Gq, and Gm, as well as 
results for the deuteron's static properties /x^, Qd, and the deuteron charge radius. We 
conclude in Section 



2 The deuteron current 

The heavy-baryon chiral perturbation theory (HB^PT) Lagrangian is organized accord- 
ing to the powers of P which appear in the classical Lagrange density. The pieces of 
the leading-order (0(P)) heavy-baryon Lagrangian relevant to the computation to be 
presented here are: 

= N\w ■ D)N + qaN^u ■ SN, (14) 



with: 
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and V chosen to be w = (1, 0, 0, 0), so that: 



2 

We also choose the pion interpolating field such that 



5 = (0,S); S = ^, (17) 
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Afj_ is the photon field. Note that we have omitted some terms that are of higher-order in 
the pion field than we need for our calculation. 

The part of relevant for our calculation is the photon-nucleon piece. There we 
focus on the vertices, suppressed by order p,Q/M, that govern the coupling of El and 
Ml photons to the nucleon j^: 

= [(^ ■Df-D.D]N- ^N^[S„ S,] [(1 + + (1 + F^^^N, (19) 

with Ffj^i, the electromagnetic field strength tensor. Kg and are the isoscalar and isovec- 
tor parts of the anomalous magnetic moment of the nucleon. These are known experi- 
mentally, and have the values —0.12 and 3.90, respectively. 

There is also an important term whose coefficient is entirely determined by reparam- 
eterization invariance. It occurs after the Foldy-Wouthysen transformation is used to 
eliminate the lower-component of the heavy-baryon field [HU] : 

^FW = -N^'^{S-D,vu}N. (20) 

Employing the definitions above, then reorganizing the result by eliminating total deriva- 
tives and using the nucleon equation of motion, leads to the piece relevant for our study: 

^f^N = ^^^^"^'^((^ ■ 9> ■ ^)^- (21) 
The first occurrence of the finite electric radius of the isoscalar nucleon occurs in chiral 

(3) 

perturbation theory as a coefficient in the Lagrangian C)y/f. Similarly, the magnetic radius 

of the nucleon appears as a coefficient in C^^}f. In both of these Lagrangians one also 
encounters terms arising from relativistic corrections to the single- nucleon four-current. 
The coefficients of these structures are determined by reparameterization invariance, and 
can be found by taking the relativistic current operator and using the standard procedure 
for generating the non-relativistic one-body current operator as an expansion in powers 
of p/M and Q/M (see, for instance jHlEIl)- 

Finally, in C^^^ we encounter a two-nucleon operator representing a magnetic photon 
coupling to the A^A^ system 



= -^eL2iN^[S,, S.]F>^-'N)iN^N). (22) 

This short-distance two-body current will modify the magnetic moment of deuterium. 

Similarly, in C\^^i^ there is an operator which represents a quadrupole (E2) photon 
coupling to the NN system j^H] and so modifies the deuteron quadrupole moment. At 
the same order there is also an operator which modifies the deuteron charge radius P^ . 

The vertices derived from the Lagrangians (|14 p - (j22p are then used to draw all possible 
Feynman diagrams contributing to the process 'y*NN —>■ NN. A particular Feynman 
diagram then leads to an operator appearing in the sum (jl3|) . The power of P that this 
operator possesses is defined by considering all parts of the amputated Feynman diagram 
representing it, and multiplying together the "P-scaling factors" of these separate pieces. 
These factors are defined as follows: 
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• A vertex from Lr^^^ contributes a factor of P". 

• A vertex involving a photon from L^^}^, L'^^Jn^ ^^^nn contributes a factor of 

pn-l 

• Each pion propagator contributes a factor of P~^. 

• Each nucleon propagator contributes a factor of P~^. 

• A two-body graph has an additional factor of P^. 

• Each loop contributes a factor of P^. 

We now discuss the charge and current operators in turn. Such a decomposition is, 
of course, not Lorentz invariant, so here we make this specification in the Breit frame, 
where the three-momentum of the deuteron and the nucleons is as shown in Fig. ^ 




-p - q/4 



Figure 1: Three momenta of the deuteron, photon, and nucleons in the Breit frame for 
a generic one-body contribution to J^. This frame is chosen because in it the photon is 
purely space-like: q = (0, q). Time runs from right to left. 

Deuteron charge The vertex from £^^}f which represents an Aq photon coupling to the 
nucleon gives the leading-order (LO) contribution to Jq: 

= |e|. (23) 

This is depicted in Fig. Efa). 

(3) 

The most important correction to Jq arises from the insertions in which generate 
the nucleon's isoscalar charge radius. This gives a result for Jq through O(eP^): 

^oLlcture = |e|(l- J (rDg^), (24) 
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where {r'^^) isoscalar charge radius of the nucleon, for which we adopt the value: 

(r|,) = (0.777 fm)l (25) 

(Note: = holds in the Breit frame.) 

Also present at this order are relativistic corrections to the single-nucleon charge opera- 
tor. To generate the "intrinsic" current operator which can be inserted between deuteron 
wave functions calculated in the two-nucleon center-of-mass frame we employ the formal- 
ism of Adam and Arenhovel, as described in Ref. The relativistic corrections then 
fall into two categories: corrections coming from the expansion of the relativistic single- 
nucleon current in powers of p/M, and corrections due to the necessity of boosting the 
deuteron wave function from the frame where P = to the frame where P = ±q/2. 






0(e) 0{eP^) 0{eP^) 

(a) (&) (c) 

Figure 2: Diagrams representing the leading contribution to the deuteron charge operator 
[(a)], the leading two-body contribution to Jo [(b)], and the dominant short- distance piece 
[(c)]. Solid circles are vertices from and the shaded circle is the vertex from C^'^^ 
The hatched square is a four- nucleon vertex from 



When the calculation is organized in this way the dominant "relativistic effect" for 
momentum transfers of order 500 MeV is a shift in the length of q. This "length contrac- 
tion" accounts for a portion of the boost of the deuteron wave function (for details see 
Refs. I31j). The net result is that whereas the leading-order form factor Gc can be 
represented as: 

''S^°'-i»l/(&*-(p-i)*(p). P«) 

with ip the deuteron wave function, at 0{eP^) the expression is: 

GcS„. = N/|^*-(p+^).Hp). (27) 

where t] = Q'^/{AMj) was defined above. Here we have not reproduced the terms which 
scale as p/M, and we have not included the terms from Eq. ((2H). The sole effect written is 
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the one arising from the boost of the deuteron wave function, ahhough all effects occurring 
at O(eP^) are included in our computation. 

This completes the discussion of mechanisms contributing at O(eP^), or next-to- leading 
order. At 0(eP'^) — next-to-next-to- leading order — the Lagrangian (fT^ generates a tree- 
level two-body graph with an isoscalar structure, as shown in Fig. Efb). This two-body 
contribution to J° was derived by Riska in Ref . , using an argument based on matching 
to relativistic Born graphs for pion electroproduction. Here it occurs in HBxPT as a 
natural consequence of the Foldy-Wouthysen transformation which generates the relevant 
term in C^'^\ Importantly, the nuclear effective theory also has the ability to organize 
the contribution of two-body contributions, such as this, relative to the contribution of 
one-body mechanisms. 

Straightforward application of the Feynman rules for the relevant pieces of the HB^PT 
Lagrangian gives the result for this piece of the deuteron current: 



pVf (q)|p) 



ai ■ qffa ■ (p - p' + q/2) 

+ (p - p' + q/2)2 



+ (1 



(28) 



where p and p' are the (Breit-frame) relative momenta of the two nucleons in the initial 
and final-state respectively [HE] . 

The short-distance two-body currents that contribute to (r^) and Qd are depicted in 
Fig. I2tc). They do not give a contribution until O(eP^). This suggests that the charge 
operator is not particularly sensitive to short-distance physics, since two-body effects of 
range 1/A are suppressed by five powers of P relative to the LO result. 



Deuteron three-current The counting for the isoscalar three-vector current J was 
already considered in detail by Park and collaborators j2Zl- J begins at 0(eP). There 
the operator derived from is: 

J(i) = |e| (p + q/2)/M + z/i5(7 x q, (29) 

where p — q/4 is the momentum of the struck nucleon, as shown in Fig. ^ and fis is the 
isoscalar magnetic moment of the nucleon, whose value we take to be fj,s = 0.88|e|/(2M). 

As in the case of Jo, there are finite-size and relativistic corrections to Eq. (pUj) which 
are suppressed by two powers of P^. Thus, in this case they enter at O(eP^), and represent 
the NLO corrections to Gm- Loop graphs of the type depicted in Fig. |21^a) also enter at 
this order. However, it can be shown that the only effect of these loops on the isoscalar 
A^A^ current is to renormalize the magnetic moment of the nucleon: their isoscalar part 
does not have any dependence (an analogous argument is given for real photons in 
Ref. 123 )• 

At 0{eP'^) [NNLO] two kinds of magnetic two-body current enter the calculation. Park 
et al. have pointed out that when magnetic photons interact with deuterium there is a 
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single-nucleon 771 contact term in C\^}^j^ ||22|- The coefficient of this portion of the chiral 
Lagrangian was fixed in Ref. j2ZI using the KSRF relation and a resonance-saturation 
hypothesis. Alternatively, this coefficient could also be fixed by comparison to data — 
at least in principle. In either case this '-fuNN vertex generates a pion-range two-body 
current with a coefficient that is undetermined a priori, as shown in Fig. El^b). 

Meanwhile, a number of authors 123 IHHl EH IHH] , have pointed out that the short- 
distance two-body operator from the Lagrangian in Eq. contributes to J at O(eP^). 
It generates a "short-range" exchange-current contribution to Gm (see Fig. Of c)). Since 
this is only suppressed by relative to the leading contributions to Gm we would expect 
Gm to be markedly more sensitive to details of the short- distance physics than Gc- Given 
the presence of two undetermined parameters at NNLO in J we will only examine the 
leading and next-to-leading order predictions of the nuclear effective theory for Gm- 




0{eP^) 0{eP^) 0{eP^) 

(a) (6) (c) 

Figure 3: An O(eP^) loop diagram which ultimately does not generate g^-dependence in 
the nucleon isoscalar form factor [(a)], and two contributions of order O(eP^) [(b) and 
(c)]. The hatched circle is a vertex from £^^\r, while the shaded square is the vertex from 



3 Deuteron wave functions 

In order to define the computation completely it remains only to specify the deuteron wave 
functions which will be used for the evaluation of the matrix elements in Eqs. ©^Q- 
Here we will use four different kinds of wave function: 

1. A "strict" chiral perturbation theory wave function, as derived in Ref. jSH]- We 
generally employ the NLO wave function, with the cutoff chosen to be A = 600 
MeV. We also use Epelbaum et al.'s NLO wave function with A = 540 MeV for 
comparison. 

2. The N^LO wave function of Ref. ||40j. In this calculation a specific choice of cut- 
off is made, which allows for better accuracy in fitting A^A^ phase shifts. Certain 
relativistic corrections to the A^A^ potential are also included. 
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3. The wave functions derived in Ref. jH] by "integrating in" the one-pion exchange 
potential (OPEP) to a given radius R. These should be regarded as very simplistic 
potential models for deuterium. They are, however, designed to produce the correct 
values for the important deuteron properties As-, Ad, and B, as well as to include 
the standard non-relativistic OPEP (with the "modern" coupling constant). 

4. The deuteron wave function obtained using the Nijm93 meson-theoretic poten- 
tial ng. 

There is an important question of consistency with the current for all of these wave 
functions. In particular, it is well known that the charge contribution ()28|) is associated 
with so-called "relativistic corrections" to the one-pion-exchange potential. If Jq and the 
terms in OPEP suppressed byp^/M^ relative to the leading behaviour are derived within 
a consistent framework then the results for deuteron form factors should be unitarily 
equivalent jlHl EH ES]. In fact, the authors of Ref. did not consider "relativistic 
corrections" to one-pion exchange. They counted M ~ and so regarded the pieces of 
OPEP of relative order p^/M^ as being down by compared with the leading piece of 
the chiral A^A^ potential. Indeed, none of the wave functions listed under numbers one, 
three, and four above include any "relativistic corrections" to one-pion exchange. Clearly 
a fully-consistent treatment of the deuteron current and A^A^ potential in xPT which 
incorporates what has been learned about unitary equivalence |^ EH |^ is necessary 
if a definitive result is to be established for electron-deuteron scattering in the nuclear 
effective theory. 

Here our goal is less ambitious. We take wave functions presently on the market and 
use the expansion for the deuteron current discussed in Sec. |21to generate results for Gc, 
Gq, and Gm- The error resulting from inconsistencies in this procedure can be assessed by 
comparing the results obtained with the wave functions of Refs. jSHlEIlES to those found 
using the "Idaho" wave function of Ref. jlU]- Of the wave functions used here, only the 
"Idaho" wave function includes the effect of relativistic corrections to one-pion exchange 
of the type associated by unitary equivalence with the two-body charge contribution (j2Hl)- 

4 Results 

Strict chiral expansion The results of the leading-order (LO), next-to- leading order 
(NLO), and next-to-next-to-leading order calculations for Gc, using the NLO xPT wave 
function of Ref. [39 with A = 600 MeV, are displayed in Fig. El Also shown there are 
data from the compilation The xPT expansion for Jq appears to be converging for 
q < 600 MeV, but it is not converging to the data. As was already observed in the NLO 
calculations of Walzl and Meifiner |29j, a strict chiral expansion of J° does a poor job of 
describing data on Gc for > 0.1 GeV^. 



9 




Figure 4: Gc as calculated in the strict xPT expansion for Jq at leading, next-to-leading, 
and next-to-next-to-leading order, plotted against |q|. The experimental data is taken 
from the extraction of Ref . ^H] '■ upward triangles represent data from the T20 measurement 
of Ref. [in], open circle [13], solid circle [H], open squares [H], downward triangles [TT] . 
star [T^, solid squares [H], solid diamonds [T^ . 

The reason for this failure can be traced to the isoscalar nucleon form factor obtained 
in xPT [46j. That form factor is: 

G^Em = l-\{rl)Q\ (30) 

and it describes electron-nucleon scattering data only up to ~ 0.1 GeV^. The inclusion 
of heavy mesons in the chiral lagrangian remedies this situation somewhat |17| , but if we 
insist on a strict xPT expansion — or even include explicit Delta degrees of freedom in the 
theory [IH] — our description of electron-deuteron data will be limited by xPT's difficulty 
in describing single-nucleon isoscalar electromagnetic structure. 

Factorization A solution to this problem is provided by the factorization of Jq. Up to 
the order to which we work the deuteron charge operator can be written as the product 
of a piece that describes the current due to structureless nucleons and a nucleon-structure 
piece: 

(pVoIp) = (|e|<5(=^)(p' -V- qm + (pVf (q)|p))G'I(Q') + O(eP^). (31) 

(Relativistic corrections are not written here, but, at this order, factorization is also valid 
for them.) 
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Figure 5: Gc and Gq (in units of fm^) calculated with nucleon structure effects included 
via factorization (at NLO and NNLO, LO is as before). The NLO xPT deuteron wave 
function with A = 600 MeV was used. Legend as in Fig. 0] 



Here we focus on the ability of nuclear effective theory to describe deuteron structure, 
and so we choose to apply the chiral expansion to the ratios of form factors: 



c 



G 



and 



(32) 



To do this we compute the ratio Jq/G^, i. e. the electric response of a deuterium nucleus 
containing structureless nucleons. Then, in order to compare with the compilation P^Bj . 
the ratios ()32|) are multiplied by the parameterization of G^ found in Ref. jlH]- The 
results obtained by this procedure are shown in Fig. This time the expansion not only 
converges, provided that Q < 700 MeV, but also reproduces data on both Gc and Gq in 
this range of Q. 

Expanding the quantities (jH^ in the effective theory sidesteps xPT's problems in 
describing isoscalar nucleon structure. We find that the chiral expansion for these ratios 
is in good agreement with data. Since these are the type of quantities which must be 
calculated in order to extract nucleon-structure information from deuteron data the results 
shown in Fig. are quite encouraging in this regard. 

Turning to the magnetic form factor, factorization also holds there, to the order to 
which we work, and so we compute the chiral expansion for the ratio J'^/G'^j. Since we 
only calculate J+ to NLO it is difficult to judge the convergence of the series, but the 
description of the data is quite good over the range Q < 500 MeV. 
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0.0 200.0 400.0 600.0 800.0 



Q (MeV) 

Figure 6: The deuteron magnetic form factor as calculated to LO with the xPT NLO wave 
function (black short-dashed line) and NLO with the xPT NLO wave function (solid black 
line) and the i? = 1.5 fm + OPEP wave function (red long-dashed line). Factorization is 
used to include nucleon structure in the NLO results. Experimental data from deuteron 
magnetic moment, open triangle [7j; the parameterization of Ref. [T7j, open squares; and 
measurements of B^Q"^): solid circles jSU], open diamonds |5T], and stars |S2]- 

Estimating the size of short-distance effects In order to judge the sensitivity of 
this observable to short- distance effects, in Fig.lHlwe also show the result for Gm obtained 
with a simple short- distance -|- OPEP wave function [41 \ This wave function and the 
xPT NLO wave function differ only at distances r -C l/mir, and so the red-dashed line's 
agreement with data to Q ~ 900 MeV should be regarded as fortuitous. From an EFT 
point of view, the difference between the red-dashed and solid curves in Fig. El is a short- 
distance effect. Such effects enter at NNLO in this observable, and so the sizable impact 
of physics at distances r ~ 1/A on Gm that is seen in Fig. IHlis not surprising. 

In contrast, short-distance contributions to Gc and Gq do not occur until O(eP^). 
As with Gm, we can estimate their impact by computing the form factors with different 
deuteron wave functions — see Fig. [7| The results for Gc and Gq are largely the same for 
Q < 600 GeV. The most noticeable difference occurs around the zero of Gc — a region 
where sensitivity to details of deuteron physics is well-established. 

Intriguingly, the band representing different assumptions about short-distance physics 
is quite narrow out to values oi Q> 800 MeV when Gq is considered. This suggests 
that the shape of Gq is not strongly affected by short- distance physics, and higher-order 
corrections to it may well be small. (A similar conclusion was reached without the use of 
nuclear effective theory in Ref. |^.) 

The curves of Figs. |3HZIare not, strictly speaking, predictions of xPT for Gc and Gq. 
In terms of the chiral expansion for these form factors a particular class of higher-order 
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200.0 



400.0 600.0 

Q (MeV) 



800.0 1000.0 



Figure 7: Results with different wave functions for Gc and Gq (in units of fm^). Solid 
black and blue lines are with wave functions from Ref. [SHI, purple dot-dashed with that 
of jlO], dotted jl2], and red and green long-dashed jH]. 

terms for electron-nucleon scattering have been resummed: the class of terms responsible 
for reproducing the "experimental" G^. Nevertheless, the results of the procedure we 
have adopted show that nuclear effective theory does a good job of describing deuteron 
structure — and especially the deuteron charge distribution — out to surprisingly high mo- 
mentum transfers. 



Deuteron static properties: As far as deuteron static properties are concerned it is 
irrelevant how nucleon structure is included in the calculation. We have computed: 



-6 



dG. 



c 



dQ^ 



(33) 



when Gc is calculated with structureless nucleons. The result for = (r^)^/^ is shown in 
TablelH together with results for /x^ and Qd- Once again, the convergence of the expansion 
is very good, with the leading-order result capturing most of the physics of these static 
properties. 

In order to assess the sensitivity of these quantities to short-distance effects we have 
computed r^, /x^, and Q^, with a variety of deuteron wave functions. The results are 
summarized in Table |21 and agreement with experimental data is very good. The ~ 0.5% 
discrepancy in is certainly consistent with the expected size of the corrections 
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Order 


ra (fm) 


fid (n.m.) 


Qd (fm2) 


LO 


1.975 


0.8591 


0.2660 


NLO 


1.984 


0.8531 


0.2641 


NNLO 


1.987 


Experiment 


0.2764 



Table 1: Deuteron static properties computed with the NLO xPT wave function (A = 600 
MeV) at LO, NLO, and NNLO. At NNLO fi^ can be exactly reproduced by adjusting the 
coefficient L2 in Eq. (j22I)- (The numerical error in each quantity is ±1 in the last significant 
figure quoted.) 

omitted here, while the ~ 1% discrepancy in fid is perhaps less than one would naively 
expect, given the that effects of relative order were not included in this computation 
of the deuteron's magnetic moment. 





Experiment 


Nijm 93 


XPT NLO 
A = 600 MeV 


NNLO Idaho 


OPEP + short 
R = 1.5 fm 


As (fm-i/2) 


0.8846(8) 


0.8842 


0.869 


0.885 


0.8845 


Ad/As 


0.0256(4) 


0.0252 


0.0248 


0.0245 


0.0253 


B (MeV) 


2.224575(9) 


Fit 


2.161 


Fit 


2.2246 


Td (fm) 


1.971(5) 


1.979 


1.987 


1.984 


1.975 


fid (n.m.) 


0.857406(1) 


0.848 


0.853 


0.847 


0.847 


Qd (fm^) 


0.2859(3) 


0.280 


0.276 


0.291 


0.280 



Table 2: Deuteron static properties at NNLO (r^ and Qd) and NLO (fid) for a range of 
deuteron wave functions. At NNLO fid can be reproduced exactly. 

On the other hand, it is apparent that Qd is much more sensitive to short- distance 
physics than either rd or fid- Its value varies by about 5% between models with the same 
pion-range, but different short-distance, physics. The counterterm that would absorb this 
sensitivity is nominally of O(eP^), which we estimate to be almost ten times smaller than 
is necessary to absorb the variation seen in Table El Whether this counterterm should 
be promoted to a lower order — as has been argued in Refs. HI] — cannot be properly 
determined until higher-order calculations of Qd are performed and a systematic study of 
its renormalization-group evolution is made. 

5 Conclusion 

Chiral perturbation theory, applied to the deuteron four-current in the fashion suggested 
by Weinberg [231231211; produces an expansion in increasing powers of small momenta 



14 



(P) for the deuteron form factors Gc, Gq, and Gm- However, this expansion fails to 
reproduce the experimental data at momentum transfers Q ~ 300 MeV [211 • The failure, 
however, lies not in xPT's description of deuteron structure, but with its difficulties in 
describing isoscalar nucleon structure. Applying a chiral expansion to the ratio of deuteron 
and nucleon form factors yields NNLO results for Gc and Gq that agree with data to 
Q ~ 700 MeV. Gc and Gq are also relatively insensitive to short- distance physics over 
this range. 

The magnetic form factor, Gm, was computed up to NLO, and turns out to be more 
sensitive to short- distance physics. This result is anticipated within the effective theory, 
since short-distance two-body currents are suppressed by three powers of P relative to 
leading in Gm, but are down by two additional powers of P in Gc and Gq. 

Deuteron static properties are also well reproduced, although shows significant 
variability when different assumptions about deuteron short-distance physics are made. 
This may be associated with the "Q^-puzzle" : the inability of modern potential models 
to reproduce the experimental value for this quantity |5^. Any possible resolution of this 
puzzle within the nuclear effective theory will require the computation of higher-order 
effects in Jq, including two-pion-exchange contributions to the deuteron four-current. 
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